Observation of a resonance X(1835) in J/ip — > 777+77 rj' 
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The decay channel J/ip — > ■yiY + Tv~n' is analyzed using a sample of 5.8 x 10 7 J ftp events collected 
with the BESII detector. A resonance, the X(1835), is observed in the 7r + 7r~r;' invariant mass 
spectrum with a statistical significance of 7.7a. A fit with a Breit-Wigner function yields a mass 
M = 1833.7 ±6.1{stat)±2.7{sy st) MeV/c 2 , a width T = 67.7 ± 20. 3{stat)± 7. 7 (syst) MeV/c 2 and a 



2 



product branching fraction B(J/4> -» ■yX) ■ B(X -> tt+tt - ?/) = (2.2 ± 0.4(staf) ± 0A(syst)) x 1CT 4 . 
The mass and width of the X(1835) are not compatible with any known meson resonance. Its 
properties are consistent with expectations for the state that produces the strong pp mass threshold 
enhancement observed in the J/ip — > *ypp process at BESII. 

PACS numbers: 12.39.Mk, 13.75.Ev, 12.40.Yx, 13.20.Gd 



An anomalous enhancement near the mass threshold in 
the pp invariant mass spectrum from J/ip — > ypp decays 
was reported by the BES II experiment Q. This en- 
hancement was fitted with a sub-threshold S'-wave Brcit- 
Wigner resonance function with a mass M = 18591-^+25 
MeV/c 2 , a width V < 30 MeV/c 2 (at the 90% C.L.) and a 
product branching fraction (BF) B(J '/tjj —> jX) ■ B(X — > 
pp) = (7.0 ± 0A(stat)±H(syst)) x 10~ 5 . This surpris- 
ing experimental observation has stimulated a number 
of theoretical speculations 0, IE 0> IE IE 01 anc ^ mo- 
tivated the subsequent experimental observation of a 
strong pA mass threshold enhancement in J/ip — > pK~ A 
decay |8j. Among various theoretical interpretations of 
the pp mass threshold enhancement, the most intriguing 
one is that of a pp bound state, sometimes called bary- 
onium 0, IE 0, which has been the subject of many 
experimental searches (lfij. 

The baryonium interpretation of the pp mass enhance- 
ment requires a new resonance with a mass around 
1.85 GeV/c 2 and it would be supported by the obser- 
vation of the resonance in other decay channels. Possible 
decay modes for a pp bound state, suggested in Ref. 0,0 , 
include 7r + 7r~?/. In this letter, we report an analysis on 
the J/ip — > 77r + 7r~?/ decay channel, and the observa- 
tion of a resonance in the 7r + 7r _ r/' mass spectrum with 
a mass around 1835 MeV/c 2 , where the rj meson is 
detected in two decay modes, rf — > ir + ir~r/(r/ — > 77) 
and 7/ — > 7/). In the following, this resonance is des- 
ignated as the X(1835). The results reported here are 
based on a sample of5.8x 10 7 J/ip decays detected with 
the upgraded Beijing Spectrometer (BESII) at the Bei- 
jing Electron-Positron Collider (BEPC). 

BESII is a large solid-angle magnetic spectrometer 
that is described in detail in Ref. [llj. Charged parti- 
cle momenta are determined with a resolution of a p /p = 

1.78% v /l+p 2 (GeV/c 2 ) in a 40-layer cylindrical main 
drift chamber (MDC). Particle identification is accom- 
plished by specific ionization (dE/dx) measurements in 
the MDC and time-of-flight (TOF) measurements in a 
barrel-like array of 48 scintillation counters. The dE/dx 
resolution is a 'dE/dx = 8.0%; the TOF resolution is mea- 
sured to be cftof = 180 ps for Bhabha events. Out- 
side of the time-of-flight counters is a 12-radiation-length 
barrel shower counter (BSC) comprised of gas tubes 
interleaved with lead sheets. The BSC measures the 
energies and directions of photons with resolutions of 
a E /E ~ 21%/ v /E(GeV), cr = 7.9 mrad, and a z = 
2.3 cm. The iron flux return of the magnet is instru- 



mented with three double layers of counters that are used 
to identify muons. In this analysis, a GEANT3-based 
Monte Carlo (MC) package with detailed consideration 
of the detector performance is used. The consistency be- 
tween data and MC has been carefully checked in many 
high-purity physics channels, and the agreement is rea- 
sonable |l2j |. 

For the J/ip — ^ 771"^ ~tt~ rj '(?/ — » 7r + 7r — r/, 17 — > 77) chan- 
nel, candidate events are required to have four charged 
tracks, each of which is well fitted to a helix that is within 
the polar angle region | cos 8\ < 0.8 and with a transverse 
momentum larger than 70 MeV/c. The total charge of 
the four tracks is required to be zero. For each track, 
the TOF and dE/dx information are combined to form 
particle identification confidence levels for the 7r, K and 
p hypotheses; the particle type with the highest confi- 
dence level is assigned to each track. At least three of 
the charged tracks are required to be identified as pions. 
Candidate photons are required to have an energy deposit 
in the BSC greater than 60 MeV and to be isolated from 
charged tracks by more than 5°; the number of photons 
are required to be three. A four-constraint (4C) energy- 
momentum conservation kinematic fit is performed to the 
7r + 7T~7r + 7r~777 hypothesis and the \\c ^ s required to be 
less than 8 and also less than the % 2 for the kinemati- 
cally similar K + if ~7r + 7r _ 777 hypothesis. An rj signal 
is evident in the 77 invariant-mass distribution of all 77 
pairings (Fig. ^a)). In order to reduce combinatorial 
backgrounds from ir° — > 77 decays, we require that the 
invariant masses of all 77 pairings are greater than 0.22 
GeV/c 2 . Candidate 77 mesons are selected by requiring 
|M 77 — m, ( | < 0.05 GeV/c 2 . The events are then sub- 
jected to a five-constraint (5C) fit where the invariant 
mass of the 77 pair associated with the 77 is constrained 
to m^, and xlc < 15 is required. The 5C-fit improves 
the M^+^-jj mass resolution from 20 MeV/c 2 (for the 
4C-fit) to 7 MeV/c 2 . Figure [2T>) shows the 7r+7r"?7 in- 
variant mass distribution after the 5C fit, where a clear 
r/ signal is visible. For rj' candidates, we select n + -K~ri 
combinations with \M 7r +^- v — m,f \ < 0.015 GeV/c 2 . In 
a small fraction of events, more than one combination 
passes the above selection. In these cases the combina- 
tion with Af^+jr-jj closest to rj' mass is used The 
7r + 7r~?/ invariant mass spectrum for the selected events 
is shown in Fig. ^c), where a peak at a mass around 
1835 MeV/c 2 is observed. 

For the J/ip — » jtt^ ~tt~ rj" '(7/ — * jp) channel, events 
with four charged tracks (with zero net charge) and two 
photons are selected. At least three of the charged tracks 
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FIG. 1: Invariant mass distributions for selected J/ip — ► 
^■k^ rf (rf — > ■K + -K^rj,r) — > 77) candidate events, (a) The 
invariant mass distribution of 77 pairs, (b) The -K + ir~r] in- 
variant mass distribution, (c) The 7t + 7t - 7/ invariant mass 
distributions; the open histogram is data and the shaded his- 
togram is J/ip — > 77r + 7r~7/ phase-space MC events (with ar- 
bitrary normalization) . 
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FIG. 2: Invariant mass distributions for the selected J/ip — ► 
7r + 7r" 77' (77' — ► 7p) candidate events, (a) The invariant mass 
distribution for 77 pairs, (b) The 77r + 7r~ invariant mass 
distribution, (c) The -K + iv~r]' invariant mass distributions: 
the open histogram is data and the shaded histogram is from 
J/ip — > r y-K + -K~rj' phase-space MC events (with arbitrary nor- 
malization) . 



are required to be identified as pions. These events are 
subjected to a 4C kinematic fit to the 7r + 7r~7r + 7r~77 hy- 
pothesis, and the xlc ^ s required to be less than 8 and less 
than the x 2 for the K + K~7r + ir~j-f hypothesis. At this 
stage of the analysis, the primary remaining background 
contributions are due to J/ip — > tt + ir~ tt + ix~ tt° , J/ip —> 
tt + tt~ tt + tt~ r\ , and J/ip — > lo(lo ~ * jtt )^ it~ ir + it~ ; these 
produce peaks at m^o, m v and to w in the 77 invari- 
ant mass distribution shown in Fig. [2{a). We suppress 
these backgrounds by rejecting events with M 77 < 0.22 
GeV/c 2 , |M 77 - m v \ <0.05 GeV/c 2 or 0.72 GeV/c 2 < 
M 77 < 0.82 GeV/c 2 . To select p and if signals, all 
7r + 7r~ and 77r + 7r _ combinations are considered. The 
77r + 7r~ invariant mass distribution shows an rf signal 
(Fig.Etb)). We require that \M„+„- -m p \ < 0.2 GeV/c 2 
and |M 77r +,r- — m n '\ < 0.025 GeV/c 2 . If more than one 
combination passes these criteria, the combination with 
M~„+ n - closest to m v t is selected [l^. For this channel 
there is also a distinct peak near 1835 MeV/c 2 in the 
tt + tt~i]' invariant mass spectrum (Fig. Pic)). 

To ensure that the peak near 1835 MeV/c 2 is not due 
to background, we have made extensive studies of poten- 
tial background processes using both data and MC. Non- 



77' processes are studied with rf mass-sideband events. 
The main background channel, J/ip — > ir°ir + Tr~r]' , and 
other background processes with multi-photons and/or 
with kaons are reconstructed with the data. In addition, 
we also checked for possible backgrounds with a MC sam- 
ple of 60 million J/ip decays generated by the LUND 
model [l4| . None of these background sources produce a 
peak around 1835 MeV/c 2 in the 7r + 7r~?/ invariant mass 
spectrum. 

Figure 13 shows the Tr + ir~r]' invariant mass spectrum 
for the combined J/ip — > 771+ n~ rf \rf — > i: + n~rf) and 
J/ip — > 77r + tt~ 'rf '{if — > 7/9) samples (i.e., the sum of the 
histograms in Figs.^c) and[2Jc)). This spectrum is fit- 
ted with a Breit-Wigner (BW) function convolved with a 
Gaussian mass resolution function (with a — 13 MeV/c 2 ) 
to represent the X(1835) signal plus a smooth polyno- 
mial background function. The mass and width obtained 
from the fit (shown in the bottom panel of Fig. 0) are 
M = 1833.7 ±6.1 MeV/c 2 and V = 67.7 ±20.3 MeV/c 2 . 
The signal yield from the fit is 264 ± 54 events with 
a confidence level 45.5% ( x 2 /tio./. = 57.6/57) and 
— 21nL = 58.4. A fit to the mass spectrum without a 
BW signal function returns — 21nL = 126.5. The change 
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in — 2\nL with A(d.o.f.) = 3 corresponds to a statistical 
significance of 7.7 a for the signal. 

Using MC-determined selection efficiencies of 3.72% 
and 4.85% for the rf — ► ir + ir~rj and rj — > 7/9 modes, 
respectively, we determine a product BF of 

£(J/V -> 7 X(1835)) • £(X(1835) -► tt+tt-j/) = 
(2.2 ± 0.4) x 10~ 4 . 
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FIG. 3: The iv + ir~r]' invariant mass distribution for selected 
events from both the J/ip — > "/-k + ir~ r)' (r]' — > -K + ir~r],r] — > 77) 
and J/-!/ 1 — > n~ r]' (ji' — > 7p) analyses. The bottom panel 
shows the fit (solid curve) to the data (points with error bars); 
the dashed curve indicates the background function. 

The consistency between the two rf decay modes is 
checked by fitting the distributions in Fig. (He) and 
Fig- I3c) separately with the method described above. 
The fit to Fig. Etc) gives M = 1827.4 ± 8.1 MeV/c 2 and 
r = 54.2 ± 34.5 MeV/c 2 with a statistical significance of 
5.1cr. From the 68±26 signal events obtained from the fit, 
the product BF is B(J/tp -> 7 X(1835)) • B(X(1835) -> 
TT + 7r~r]') = (1.8±0.7)xl0 -4 . Similar results are obtained 
if we only apply 4C kinematic fit in this analysis. For the 
fit to Fig. |2{c) , the mass and width are determined to be 
M = 1836.3 ± 7.9 MeV/c 2 and T = 70.3 ± 23.1 MeV/c 2 
with a statistical significance of 6.0 a. For this mode 
alone, the signal yield of 193±43 signal events correspond 
to B(J/if> -> 7 A(1835)) ■ £(X(1835) -> tt+tt-?/) = 
(2.3 ± 0.5) x 10~ 4 . The X(1835) mass, width and prod- 
uct BF values determined from the two rf decay modes 
separately are in good agreement with each other. 

The systematic uncertainties on the mass and width 
are determined by varying the functional form used to 
represent the background, the fitting range of the mass 
spectrum, the mass calibration, and possible biases due 
to the fitting procedure. The latter arc estimated from 
differences between the input and output mass and width 
values from MC studies. The total systematic errors on 



the mass and width are 2.7 MeV/c 2 and 7.7 MeV/c 2 , 
respectively. The systematic error on the branching frac- 
tion measurement mainly comes from the uncertainties 
of MDC simulation (including systematic uncertainties of 
the tracking efficiency and the kinematic fits), the photon 
detection efficiency, the particle identification efficiency, 
the rf decay branching fractions to ir + Tr~i] and 7 p, the 
background function parameterization, the fitting range 
of the mass spectrum, the requirements on numbers of 
photons, the invariant mass distributions of 77 pairs in 
the two analyses, the 7r + 7r~ invariant mass distribution 
in rf — > 77T+7T - decays, MC statistics, the total num- 
ber of J /if) events [l5j, and the unknown spin-parity of 
the X(1835). For the latter we use the difference be- 
tween phase-space and a J PC = h hypothesis for the 
X(1835). The total relative systematic error on the prod- 
uct branching fraction is 20.2%. 

In summary, the decay channel J /if) — > 7 7r + 7r~?/ is 
analyzed using two rf decay modes, rf — ► ■k + tx~t} and 
rf — > 7 p. A resonance, the X(1835), is observed with 
a high statistical significance of 7.7er in the 7r + 7r _ r/ 
invariant mass spectrum. From a fit with a Breit- 
Wigner function, the mass is determined to be M = 
1833.7 ± 6.1(stat) ± 2.7{syst) MeV/c 2 , the width is T = 
67.7 ± 20.3(sia<) ± 7.7{syst) MeV/c 2 , and the product 
branching fraction is B(J/if> — ► jX) ■ B(X —> 7T + 7r~?/) 
= (2.2 ± OA(stat) ± OA(syst)) x 10~ 4 . The mass and 
width of the X(1835) are not compatible with any known 
meson resonance [l6j. In Ref. [16|, the candidate clos- 
est in mass to the X(1835) is the (unconfirmed) 2 ^ 
ry 2 (1870) with M = 1842±8 MeV/c 2 . This state's width, 
T = 225 ± 14 MeV/c 2 , is considerably larger than that 
of the X(1835) (see also IT], where the 2 + compo- 
nent in the r/irir mode of J / if) radiative decay has a mass 
1840 ± 15 MeV/c 2 and a width 170 ± 40 MeV/c 2 ). 

We examined the possibility that the X(1835) is re- 
sponsible for the pp mass threshold enhancement ob- 
served in radiative J /if) — > 7 pp decays 0. It has been 
pointed out that the 5-wave BW function used for the 
fit in Ref. Q should be modified to include the effect of 
final-state-interactions (FSI) on the shape of the pp mass 
spectrum 0, 0] • Redoing the S'-wave BW fit to the pp 
invariant mass spectrum of Ref. E| including the zero 
Isospin, S'-wave FSI factor of Ref. |7| , yields a mass M — 
1831 ± 7 MeV/c 2 and a width T < 153 MeV/c 2 (at the 
90% C.L.) |l8| ; these values are in good agreement with 
the mass and width of X(1835) reported here. Moreover, 
according to Ref. , the iririf decay mode is expected to 
be strong for a pp bound state. Thus, the X(1835) reso- 
nance is a prime candidate for the source of the pp mass 
threshold enhancement in J/tp — > ^ypp process. In this 
case, the J PC and I G of the X(1835) could only be 0"+ 
and + , which can be tested in future experiments. Also 
in this context, the relative pp decay strength is quite 
strong: B{X -> pp)/B{X -> tt+tt-?/) - 1/3 19]. Since 
decays to pp are kinematically allowed only for a small 
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portion of the high-mass tail of the resonance and have 
very limited phase space, the large pp branching fraction 
implies an unusually strong coupling to pp, as expected 
for a pp bound state [2(|- However, other possible inter- 
pretations of the A(1835) that have no relation to the pp 
mass threshold enhancement are not excluded. 
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